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Broadband modeling of GRB 021004 


Jeremy S. Heyl^’^ & Rosalba Perna^’^ 

ABSTRACT 

We present a broadband modeling of the afterglow of GRB 021004. The 
optical transient of this burst has been detected very early and densely sampled 
in several bands. Its light curve shows signihcant deviations from a simple power 
law. We use the data from the X-ray to the /—band gathered in the hrst month 
of observations, and examine three models. Two involve variations in the energy 
of the shock. The hrst (energy injection) allows only increases to the shock 
energy, while the second (patchy shell) allows the energy to increase or decrease. 
In the hnal model (clumpy medium) the energy of the shock is constant while 
the density varies. While all three models reproduce well the optical bands, the 
variable density model can best account for the X-ray data, and the energy- 
injection model has the poorest ht. None of the models can account for the 
modest color variations observed during the hrst few days of the burst. 

Subject headings: gamma rays: bursts 


1. Introduction 

The gamma-ray burst GRB 021004 was detected by HETE II at 12:06 UT on the 4th 
of October 2002 (Shirasaki et ah 2002). Observations after about 9 minutes from the trigger 
revealed a fading optical transient (Fox 2002), which was densely sampled in several bands, 
especially at early times. An X-ray counterpart has also been reported by Sako & Harrison 
( 2002 ). 

The afterglow of GRB 021004 has shown several unusual features. Absorption lines from 
high velocity material have been detected in the spectrum (Mirabal et ah 2002; Salamanca 
et ah 2002; Moller et ah 2002), evidence of short time scale variability has been presented 
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(Bersier et al. 2002), and the light curves in various bands display signihcant deviations from 
a simple power law decay, with several bumps and flattenings. A hrst, pronounced bump 
around 10^ — 10^ sec was modeled by Lazzati et al. (2002) with a density contrast in the 
ambient medium density, while Kobayashi & Zhang (2002) interpreted it as due to emission 
from the reverse shock. Nakar et al. (2002) htted the R band data and found that several 
models provide good fits to the data, both with varying density (at fixed energy) and with 
varying energy (at hxed density). Determining the appropriate physical model is crucial in 
order to understand what the progenitors of GRBs are, and the large wealth of data available 
for GRB 021004 in several bands allows detailed modeling. 

In this Letter, we present a modeling of the light curve of GRB 021004 from the X- 
ray to the /—band. We hnd that a model in which the bumps are produced by density 
fluctuations (and not by energy injection) best accounts for both the optical and the X-ray 
data. We derive the density prohle that best hts the data and discuss its implications for 
GRB progenitor scenarios. 


2. Models 

To calculate the emission from the afterglow we have used the models of Sari et al. (1998) 
as extended by Nakar et al. (2002) for a varying density or energy"^. Although both the energy 
of the afterglow and the density of the surrounding medium can vary, we will consider a GRB 
shock whose energy changes while propagating through a wind with a constant mass-loss rate 
or a const ant-density ISM. The calculations show that the former model requires the energy 
of the shock to increase monotonically so we will call it the “the energy-injection model”. 
The latter cases requires a shock whose energy varies up and down so we will call it “the 
patchy-shell model”. The third model we consider is a shock of constant energy traveling 
through medium with varying density (the clumpy-medium model). 

We refer the reader to Nakar et al. (2002) for the necessary equations, and we summarize 
their results here. If < Vm (fast cooling), the flux density is proportional to between 
Vc and Vrni to below Vc and to above Vm- In both models for the assumed values of 
the electron and magnetic held energy fractions of 0.1 and 0.01, fast cooling lasts for about 
the hrst thousand seconds after the burst; therefore, only the data of Fox (2002) lies in the 
fast-cooling regime. If Vra < (slow cooling), the hux density goes as v~L>-Lt‘^ between 


"‘Note that these models assume that all of the emission originates along the line of sight. The emission 
actually comes from a surface of equal arrival times, so the actual lightcurve should be somewhat smoother 
than predicted by these models. 
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and Uc, below and above Vc- Vm passes through the optical and near infrared 

a few thousand seconds after the burst. On the other hand, changing the energy of the 
visible portion of the afterglow either through energy injection or a patchiness changes the 
flux density proportionally at all frequencies. 


3. Results 

We compare the three models with the broadband optical culled from the GCN and the 
literature for the first day after the burst (Fox 2002; Uemura et ah 2002; Weidinger et al. 
2002; Oksanen et al. 2002; Matsumoto et al. 2002; Barsukova et al. 2002; Klotz et al. 2002; 
Holland et al. 2002; Pandey et al. 2002), the X-ray data obtained by Sako & Harrison (2002) 
and reanalyzed ourselves, and the later optical data from Bersier et al. (2002), Holland et al. 
(2002) and Pandey et al. (2002) . We adjust either the energy injection or the density of the 
surrounding medium to £t the R-band data. The Chandra observations of Sako & Harrison 
(2002) constrain the electron-injection power-law exponent p to be about 2.2. Matheson 
et al. ( 2002 ) give a redshift of 2.323 for the optical transient. Assuming the cosmographic 
parameters VLm = 0.3, Ca = 0.7 and Hq = 70 km s“^Mpc“^ yields a luminosity distance of 
5.8 X 10^® cm to the gamma-ray burst. This distance results in an isotropic energy for the 
afterglow of about 9 x 10^^ ergs during the Chandra observation Sako & Harrison (2002). 
We account for the cosmological time dilation and A:—correction. 

Because we varied either the energy of the forward shock or the density of the sur¬ 
rounding medium to obtain a good fit to the observations in the R—band and adjusted the 
extinction to fix the average flux in the other bands, the key diagnostic between the two 
models is that density variations only affect the flux between the cooling frequency (Lazzati 
et al. 2002; Nakar et al. 2002) and Vm and are achromatic in this regime. Variations in 
the energy change the flux achromatically throughout the near infrared, optical and X-ray 
regimes. 

Fig. 1 compares the two energy-variation models with the data from the first month 
after the burst. For the energy-injection model we have assumed a wind density profile 
that corresponds to a mass-loss rate of 1.3 x 10“® M 0 yr“^, assuming a wind velocity of 
1000 km s“^. The local extinction is Ay/ = 0.34, where V is the local V—band. The patchy- 
shell model has a constant density medium with n = 0.1 cm“^ and a local extinction of 
Ayi = 0.33. Fig. 2 depicts the clumpy-medium model. In this model, the isotropic energy of 
the ejecta is 8.9 x 10^^ ergs and the local extinction is Ay/ = 0.33. This energy agrees with 
the isotropic energy in gamma rays of the CRB itself (Lamb et al. 2002), 2.06 x 10^^ ergs if 
one assumes that twenty-percent of the energy in the ejecta is converted to prompt gamma 
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Fig. 1.— The panels show the results of the two energy-variation models superimposed on 
the observations. The left panel depicts the energy-injection model: the energy only increases 
as the shock propagates and the density drops as r“^. The right panel shows the patchy-shell 
model: the density is constant and the energy varies up and down. From bottom to top, 
the curves trace the X-ray flux at 2 keV, U + 2, B + 1, V + 0.5, R, and I — 0.5. The X-ray 
magnitude is dehned using a zero point 800 Jy at 2 keV. 
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Fig. 2.— The panel shows the results of the clumpy-medium model. Here the energy of the 
shock is constant but the density of the medium is allowed to vary 
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rays (Frail et al. 2001). For all three models the Galactic value of Ay is taken to be 0.21 in 
the direction of the burst (Bersier et al. 2002). We used the Milky-Way and LMC extinction 
models of Pei (1992) to model the Galactic and host dust respectively. For comparison 
Holland et al. (2002) hnd Ay^ = 0.26 ± 0.04 using Pei’s SMG model. 

All three models provide excellent hts to the optical and near infrared data from one 
hour until about six days after the burst. The patchy-shell model fails to predict the flux 
during the first hour of the burst and overpredicts the flux from the afterglow after six days. 
The initial points may be explained by emission from a reverse shock (e.g. Pandey et al. 
2002) or an increase in the density of the medium near the progenitor. Holland et al. (2002) 
argue that the jet break occurs at about six days, so because our models assume an isotropic 
afterglow, the presence of a jet can account for this discrepancy. The energy-injection and 
clumpy-medium models can £t the entire i?—band light curve up to one month after the 
burst. 

Fig. 3 depicts the required variation in the energy of the visible portion of the fireball 
in the energy-variation models and Fig. 4 depicts the run of density in the clumpy-medium 
model. In the clumpy-medium and patchy-shell models, the required variations are modest 
so we are justified in neglecting the contribution of reverse shocks to the emission (Lazzati 
et al. 2002). The energy-injection model requires the energy of the shock to increase by a 
factor of forty, so the simple models of Nakar et al. (2002) may not be applicable. Note 
that, in the clumpy-medium model, the rise of the density at small radii derives from fitting 
the early optical data (Fox 2002) such a rise would not be required if that emission were 
originated in the reverse shock as claimed by (Pandey et al. 2002) or if the early data were 
also in the slow cooling regime due to a different choice of the parameters Ce and €b (Lazzati 
priv. comm). Although the transition from fast cooling to slow cooling is dependent on Ce 
and es, during either regime changes in and es are degenerate with changes in the energy 
of the shock and the density of the medium. 

The largest difference between the models is at frequencies above Uc- The Ghandra 
observations (Sako & Harrison 2002) are contemporaneous with optical and near infrared 
measurements (Bersier et al. 2002; Holland et al. 2002; Pandey et al. 2002). The optical/NIR 
observations show significant variation from power-law decay during this epoch. The X- 
ray observations are consistent with the power-law decay In the clumpy-medium 

model, the X-ray flux follows a strict power-law decay during this epoch because the cooling 
frequency lies below the X-rays. On the other hand, the energy-variation model results in 
complicated behavior in the X-rays which is absent in the data. 

We have reanalyzed the data obtained by Sako & Harrison (2002) to obtain better 
constraints on the departures to the power-law evolution of the flux in the X-ray band 
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Fig. 3.— The isotropic energy of the afterglow as a function of radius for the energy-variation 
models. The short-dashed curve follows the energy-injection model, and the long-dashed 
curve follows the patchy-shell model. The horizontal lines show the range of radii covered 
by the observations. 
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Fig. 5.— The figure shows the value of the bulk Lorentz factor as a function of the time after 
the burst. The solid curve traces the solution with clumpy-medium model, the short-dashed 
lines follow the energy-injection result, and the long-dashed line gives the patchy-shell result. 
The horizontal lines show the range of radii probed by the observations in each of the models. 
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as they are crucial for discriminating among the various models. We looked at only the 
zeroth order image which Sako & Harrison (2002) have argued is consistent with the HETG 
spectrum. We £t the ACTS spectrum for the entire observation with an absorbed power- 
law energy distribution and obtained values consistent with the results of Sako & Harrison 
(2002). We then compared the distribution of photon arrival times for all those photons 
above the median energy with those below the median energy. Both these distributions 
are consistent with the distribution of all of the photon arrival times. This indicates that 
there was little spectral evolution during the observation, so we divided the observation into 
four intervals each with the same number of counts and £t the resulting spectrum with 
absorbed power-law energy distributions with the power-law index and absorption column 
hxed, yielding the points depicted in Figs. 1 and 2. 

To test our models robustly we used them to predict the distribution of photon arrival 
times during the X-ray observations. We compared these distributions with the observed 
distribution using the Kolmogorov-Smirnov statistic. We found that the density-variation 
model is preferred (K-S probability of 3%) over the patchy-shell model (K-S probability of 
0.3%). The energy-injection model was highly disfavored with a K-S probability of 2 x 10“^. 

If we examine in further detail, the color of the emission in the optical and near infrared 
varies significantly from day to day (Bersier et al. 2002; Matheson et ah 2002). Because all of 
these bands he below the cooling frequency, both the density-variation and energy-variation 
models cannot reproduce these color changes. This is apparent from Figs. 1 and 2. 


4. Discussion 

We find that the broadband emission from GRB 021004 can be well described by a 
relativistic fireball propagating into a clumpy medium. This model can account for the 
observed variability in the optical and near infrared while producing a power-law decay in 
the X-rays. Varying the energy in the observed portion of the fireball results in a variation 
in the X-ray flux as well which was not seen. 

The patchy-shell model provides a poorer explanation. In this model, the mean energy 
of the observed portion of the fireball can both increase or decrease (Nakar et al. 2002) as 
it travels through a constant density medium. This model can account for the near infrared 
through ultraviolet data, but it predicts a departure from power-law decay in the X-rays. 
It cannot account for the earliest data or for the fluxes more than five days after the burst. 
This time coincides with the jet break found by Holland et al. (2002). Because our models 
do not account for jet emission, the jet provides a natural explanation for the discrepancy. 



In both the clumpy-medium and the patchy-shell models F 10 at this time. Using the 
formula of Frail et ah (2001) gives a jet-opening angle of 0.09 consistent with the value of 
r. This yields a total energy of the ejecta of approximately 4 x 10^'^ erg and 9 x 10^® erg in 
the gamma-ray burst itself, somewhat less than the mean found by Frail et al. (2001). 

The least likely model is energy injection. Because of the steep decline in the second 
half of the first day, a wind density model is required to £t the data in this scenario. Without 
some significant energy injection during this period, the model underpredicts the flux later 
on. Furthermore, the energy-injection model exhibits dramatic evolution in the X-rays as 
well. This was not seen. 

We find that the afterglow of GRB 021004 was most likely produced in a clumpy wind 
flowing into a low-density, typical ISM medium. The mean density during most of the 
afterglow phase is about 0.1 cm“^. If this value is typical also at scales larger larger than 
those in which the afterglow is produced, no variability in the absorption lines in the spectra 
(associated with the host galaxy) should be expected (Perna & Loeb 1998). The typical 
density that we infer could be very well typical of mergers of two compact objects as found 
in numerical simulations (Perna & Belczynski 2002), but in this case we do not expect a 
particularly bumpy environment (besides the small scale fluctuations due to ISM turbulence, 
see Wang & Loeb 2000), while this seems to be the case for GRB 021004. Therefore, the 
association with the collapse of a massive star is a more likely scenario. A lumpy, low-density 
background could suggest an explosion which occurred inside the remnant of a previous 
supernova explosion, such as in the Supranova model of (Vietri & Stella 1998). 

Besides the magnitude of the density, the density prohle itself in the immediate envi¬ 
ronment of the burst can be be used to discern among various progenitor models (see e.g. 
Lazzati et al. 2001). We hnd that, if we fit the early optical data points within the same 
forward shock model that is used to model the later emission®, then the overall density profile 
resembles that of a wind from a massive star, as one would expect in the collapsar scenario 
MacFadyen et al. (2001). We have computed the mass loss rate that is needed to produce 
the derived density assuming a wind velocity of 1000 km s“^. Although the density that we 
hnd is somewhat lower than that assumed by Ramirez-Ruiz et al. (2001), it is typical for 
high-mass stars at subsolar metallicity (Kudritzki 2002). 

GRB 021004 exhibited several color changes during its decay (Bersier et al. 2002; Math- 
eson et al. 2002). None of the models that we have considered can explain these changes. 
Furthermore, it is impossible to reproduce the color change within the broken-power-law 
spectral model of Sari et al. (1998) even by integrating over an heterogeneous fireball {e.g. 


®This might not necessarily be the case according to (Pandey et al. 2002). 
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Nakar et al. 2002), unless the inhomogeneities are sufficiently large to shift the values of Vc 
or Vm across the bands of interest. Because the latter possibility is unlikely, we conclude that 
they are most likely either external to the hreball or a result of the approximations used in 
current afterglow models. 
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